Dynamic Stability of the 3D Axi-symmetric 
Navier-Stokes Equations with Swirl 

Thomas Y. Hou* Congming Li ^ 

February 2, 2008 

Abstract 

In this paper, we study the dynamic stabihty of the 3D axisymmetric Navier- 
Stokes Equations with swirl. To this purpose, we propose a new one-dimensional 
(ID) model which approximates the Navier-Stokes equations along the symmetry axis. 
An important property of this ID model is that one can construct from its solutions a 
family of exact solutions of the 3D Navier-Stokes equations. The nonlinear structure 
of the ID model has some very interesting properties. On one hand, it can lead to 
tremendous dynamic growth of the solution within a short time. On the other hand, it 
has a surprising dynamic depletion mechanism that prevents the solution from blowing 
up in finite time. By exploiting this special nonlinear structure, we prove the global 
regularity of the 3D Navier-Stokes equations for a family of initial data, whose solutions 
can lead to large dynamic growth, but yet have global smooth solutions. 

1 Introduction. 

Despite a great deal of effort by many mathematicians and physicists, the question of whether 
the solution of the 3D Navier-Stokes equations can develop a finite time singularity from 
a smooth initial condition with finite energy remains one of the most outstanding open 
problems ^2]- A main difficulty in obtaining the global regularity of the 3D Navier-Stokes 
equations is due to the presence of the vortex stretching, which is absent for the 2D problem. 
Under suitable smallness assumption on the initial condition, global existence and regularity 
results have been obtained for some time [TTl |HJ 120] • But these methods based on energy 
estimates do not generalize to the 3D Navier-Stokes with large data. Energy estimates seem 
to be too crude to give a definite answer to whether diffusion is strong enough to control the 
nonlinear growth due to vortex stretching. A more refined analysis which takes into account 
the special nature of the nonlinearities and their local interactions seems to be needed. 
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In this paper, we study the dynamic stabihty property of the 3D axisymmetric Navier- 
Stokes Equations with swirl. We show that there is a very subtle dynamic depletion mech- 
anism of vortex stretching in the 3D Navier-Stokes equations. On one hand, the nonlinear 
vortex stretching term is responsible for producing a large dynamic growth in vorticity in 
early times. On the other hand, the special structure of the nonlinearity can also lead to 
dynamic depletion and cancellation of vortex stretching, thus avoiding the finite time blowup 
of the Navier-Stokes equations. 

This subtle nonlinear stability property can be best illustrated by a new ID model which 
we introduce in this paper. This ID model approximates the 3D axisymmetric Navier-Stokes 
equations along the symmetry axis. By the well-known Caffarelli-Kohn-Nirenberg theory jSj 
(see also [IH]), the singularity set of any suitable weak solution of the 3D Navier-Stokes 
equations has one- dimensional Hausdorff measure zero. In the case of axisymmetric 3D 
Navier-Stokes equations with swirl, if there is any singularity, it must be along the symmetry 
axis. Thus it makes sense to focus our effort to understand the possible singular behavior of 
the 3D Navier-Stokes equations near the symmetry axis at r = 0. By expanding the angular 
velocity (u^), the angular vorticity (cu^), and the angular stream function {ip^) around r = 0, 
we obtain the following coupled nonlinear partial differential equations (see Section 2 for 
detailed derivations): 

{ul)^ + 2^1 (mi)^ = u{ui)^^ + 2 (^i)^ Ui (1) 
(cui)j + 2^i(cji)^ = z/(cui),, + (n^)^ (2) 

-{'ipl)zz = UJl, (3) 
where Ui{z,t) ^ {u^)r\r=0, ^l{z,t) ^ {uj^)r\r=0, and l/ji{z,t) ^ (lp^)j.\r=Q. 

What we find most surprising is that one can construct a family of exact solutions from 
the above ID model. Specifically, if (mi, tui, ■j/'i) is a solution of the ID model ((H)-©, then 

u^{r, z,t) = rui{z,t), iu'^{r, z,t) = ruJi{z,t), ip^ {r^ z,t) = ril)i{z,t)^ 

is an exact solution of the 3D axisymmetric Navier-Stokes equations. Thus the ID model 
captures some essential nonlinear features of the 3D Navier-Stokes equations. Further, if we 
let £t = ui, -y = —{ipi)zi and = tpi, then the ID model can be rewritten as 

{u)t + 2ij{u), = v{u),, - 2vu (4) 
{v)t + 2^({}), = v{v),, + {uf - {vf + c(t), (5) 

where ipz = —v and c{t) is an integration constant to ensure that / vdz = 0. We will 
show that if the initial value of u is small, but v is large and negative, then the solution of 
the ID model can experience large growth. On the other hand, we also find a surprising 
dynamic depletion mechanism of nonlinearities that prevents the solution from blowing up 
in finite time. This subtle nonlinear cancellation is partly due to the special nature of the 
nonlinearities, i.e. —2uv in and — in (jSJ. If one modifies the sign of the nonlinear 
term from —2uv to 2uv or changes v? — v'^ to + v'^ or even modifies the coefficient from 
—2uv to —0.9m?), the dynamic depletion mechanism can be changed completely. Another 
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interesting fact is that the convection term also helps to stabilize the solution. It cancels 
some of the destabilizing terms from the right hand side when we estimate the solution in 
a high order norm. Specifically, we find that there is a miraculous cancellation of nonlinear 
terms in the equation that governs the nonlinear quantity, u1 + v"^, i.e. 

{ul + vl)^ + 2^ {ul + vj)^ = u{ul + vl)^^ - 2u [{u^^y + iv^^f] . (6) 

Therefore, + has a maximum principle. This pointwise a priori estimate plays an 
essential role in obtaining global regularity of the ID model with or without viscosity. If 
one attempts to prove global regularity of the ID model using energy estimates, one cannot 
take full advantage of this local cancellation of nonlinearities and would run into similar 
difficulties that we encounter for the 3D Navier-Stokes equations. 

Finally, we construct a family of globally smooth solutions of the 3D Navier-Stokes equa- 
tions with large initial data of finite energy by using the solution of the ID model. Specifically, 
we look for the solution of the form: 

= r{ui{z,t)(f){r)+ui{r,z,t)) (7) 
= r{uiiz,t)(l){r)+uJiir,z,t)) (8) 
= r(7/;i(z,t)0(r) + V^i(r,z,t)), (9) 

where Ui, Ui and ipi are solutions of the ID model, 0(r) is a cut-off function to ensure 
that the solution has finite energy. By using the a priori estimate of the solution of the 
ID model and using a delicate analysis, we prove that there exists a family of globally 
smooth functions Mi(r, z, t), co'i(r, z, t) and ipii^r, z, t), such that u^, oj^ and ip^ are solutions of 
the 3D axisymmetric Navier-Stokes equations. Unlike the other known global solutions with 
small data, the solutions that we construct above using the ID model can have large dynamic 
growth for early times, which is induced by the dynamic growth of the corresponding solution 
of the ID model, but yet the solution remains smooth for all times. 

There has been some interesting development in the study of the 3D incompressible 
Navier-Stokes equations and related models. In particular, by exploiting the special structure 
of the governing equations, Cao and Titi jl] prove the global well-posedness of the 3D viscous 
primitive equations which model large scale ocean and atmosphere dynamics. By taking 
advantage of the limiting property of some rapidly oscillating operators and using non-linear 
averaging, Babin, Mahalov and Nicolaenko PP prove existence on infinite time intervals of 
regular solutions to the 3D Navier-Stokes equations for some initial data characterized by 
uniformly large vorticity. Some interesting progress has been made on the regularity of the 
axisymmetric solutions of the Navier-Stokes equations, see e.g. [U] and the references cited 
there. The 2D Boussinesq equations are closely related to the 3D axisymmetric Navier- 
Stokes equations with swirl (away from the symmetry axis). Recently, Chae ^ and Hou-Li 
have proved independently the global existence of the 2D viscous Boussinesq equations 
with viscosity entering only in the fiuid equation, but the density equation remains inviscid. 
Recent studies by Constantin-Fefferman-Majda and Deng-Hou-Yu pniE] show that the 
local geometric regularity of the unit vorticity vector can play an important role in depleting 
vortex stretching dynamically. Motivated by these theoretical results, Hou and R. Li ^3] 
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have recently re-investigated the well-known computations by Kerr JH] for two anti-parallel 
vortex tubes, in which a finite time singularity of the 3D incompressible Euler equations was 
reported. The results of Hou and Li show that there is tremendous dynamic cancellation in 
the vortex stretching term due to local geometric regularity of the vortex lines. Moreover, 
they show that the vorticity does not grow faster than double exponential in time and the 
velocity field remains bounded up to T = 19, beyond the singularity time alleged in jl6j . 
Finally, we would like to mention the recent work of Gibbon et al (see jj^] and the references 
therein) where they reveal some interesting geometric properties of the Euler equations in 
quaternion- frames. 

The rest of the paper is organized as follows. In Section 2, we will derive the ID model 
for the 3D axisymmetric Navier-Stokes equations. We discuss some of the properties of the 
ID model in Section 3 and prove the global existence of the inviscid ID model using the 
Lagrangian coordinate. Section 4 is devoted to prove the global regularity of the full ID 
model in the Eulerian coordinate. Finally in Section 5, we use the solutions of the ID model 
to construct a family of solutions of the 3D Navier-Stokes equations and prove that they 
remain smooth for all times. 

2 Derivation of the ID Model 

Consider the 3D axi-symmetric incompressible Navier-Stokes equations with swirl. 

ut + {u ■ V)u = —Vp + uAu, 

V-M = 0, (10) 

u\t=o = Uo{x), x = {x,y,z). 

Let 

= (-,^,0), ee=(--,-,o), = (0,0,1), 
\r r / \ r r / 

be three unit vectors along the radial, the angular, and the z directions respectively, r = 
a/x^ + y^. We will decompose the velocity field as follows: 

u = v^{r, z, t)er + u^{r, z, t)ee + v^{r, z, t)ez. (11) 

In the above expression, is called the swirl component of the velocity field u. The vorticity 
field can be expressed similarly 

oo = -{u^)z{r, z, t)er + uj^{r, z, t)ee + -{ru^)r{r, z, t)ez, (12) 

r 

where = vl — v^. 

To simplify our notation, we will use u and u to denote the angular velocity and vorticity 
components respectively, dropping the 6 superscript in the rest of the paper. One can derive 
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evolution equations for u and u as follows (see e.g. PU} IH]). 



Ut + V^Ur + V^Uz = Z/ ( r ] U U, (13) 



- fv^ - ^1 ^ = o;, (15) 



where t/) is the angular component of the stream function, f ^' and can be expressed in 
terms of the angular stream function tl) as follows: 

and is defined as 

V2 = 5,2 + V + 9^ (17) 

Note that equations ()13|) - (fT^ completely determine the evolution of the 3D axisymmetric 
Navier-Stokes equations once the initial condition is given. 

Now, we will derive the ID model for the 3D axisymmetric Navier-Stokes equations. By 
the well-known Caffarelli-Kohn-Nirenberg theory [Hj, the singularity set of any suitable weak 
solution of the 3D Navier-Stokes equations has one-dimensional Hausdorff measure zero. 
Thus, in the case of axisymmetric 3D Navier-Stokes equations with swirl, if there is any 
singularity, it must be along the symmetry axis, i.e. the ^-axis. Therefore, we should focus 
our effort to understand the possible singular behavior of the 3D Navier-Stokes equations 
near the symmetry axis at r = 0. 

As observed by Liu and Wang in ^H], any smooth solution of the 3D axisymmetric 
Navier-Stokes equations must satisfy the following compatibility condition at r = 0: 

m(0, z, t) = uj{0, z, t) = tp{0, z, t) = 0. (18) 

Moreover, all the even order derivatives of u, u and with respect to r at r = must vanish. 
Therefore, we expand the solution u, u and around r = as follows: 

u{r, z, t) = rui{z, t) + ^"3(2, t) + ^"5(^, t) H , (19) 

ly'-' ly''' 

Lj{r, z, t) = ruJi{z, t) + ^^3(2;, t) + ^^5(2;, t) H , (20) 

3 5 

V'(r, z, t) = r^Piiz, t) + ^^3(2;, t) + ^^5(2:, t) + ■ ■ ■ . (21) 
Substituting the above expansions into (fT^ - (fT3j) . we obtain to the leading order the 
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following system of equations: 

By canceling r from both sides and neglecting the higher order terms in r, we obtain 

^V'a + (V'l)^^)^ = ^1- 

Note that = Mrrr-(0, ^, i), = 'Wrzz(0, 2;, t). If we further make the assumption that 

the second partial derivative of iti, lo\^ ipi with respect to z is much larger than the second 
partial derivative of these functions with respect to r, then we can ignore the coupling in the 
Laplacian operator to us, 003 and -03 in the above equations. Thus, we obtain our ID model 
as follows: 

iu^)^ + (Ml), = + 2 (V^i)^ u^, (22) 

(cui), + 2iji (uji)^ = v{uj^),, + (m^)^ , (23) 

-{i^i)zz = ^1- (24) 

We remark that the above assumption implies that the solution has an anisotropic scaling, 
i.e. the solution is more singular along the ^-direction than along the r-direction. A possible 
scenario is that the solution has a pancake like structure perpendicular to the 2;-axis. 

Let u = V = — ('i/'i)^, = a;i, and = ■i/'i. By integrating the a;i equation with 
respect to z and using the relationship —-^^i/ji — coi, we can obtain an evolution equation 
for V. Now the complete set of evolution equations for u, v, and a; are given by 

{u)^ + 2i^{u), = v{u),, - 2vu, (25) 
(c^;), + 27/;(a'), = i/(6u),,+ (S^)^, (26) 

{v)^ + 2i,{v), = u{v),, + {uf - {vf + c(t), (27) 

-{ij)zz = Co, (28) 

where the constant c{t) is an integration constant which is determined by enforcing the mean 
of V equal to zero. For example, if ip is periodic with period 1 in z, then c{t) is given by 



( / v^dz — u 
Jo Jo 



c{t) = 3 / v^dz - / u^dz. (29) 
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Note that the equation for uj is equivalent to that for v. So it is sufficient to consider the 
coupled system for -u, v: 

(■u)j + 2^{u)z = l'{u)zz - S-UM (30) 
{v), + 2^{v)z = v{v)zz + {uf - {vf + c(t), (31) 

where is related to -u by -u = By (OHj) . we have -u^ = uj. 

A surprising result is that one can use the above ID model to construct a family of exact 
solutions for the 3D axisymmetric Navier-Stokes equations. This is described by the following 
theorem, which can be verified directly by substituting ()32|1 into the 3D axisymmetric Navier- 
Stokes equations and using the model equation (j22I)-(I21j)- 

Theorem 1. Let ui, ipi and ui be the solution of the ID model \2^] -\24^ and define 

u{r, z,t) = rui{z,t), uj{r, z,t) = rLJi{z,t), ip{r, z^t) = ripi{z,t). (32) 
Then {u{r,z,t), u){r,z,t), tlj{r,z,t)) is an exact solution of the 3D Navier-Stokes equations. 

Theorem 1 tells us that the ID model (j22I)-(I13) preserves some essential nonlinear struc- 
ture of the original 3D axisymmetric Navier-Stokes equations. As we will see later, the 
nonlinear structure of the ID model plays a critical role in stabilizing the solution for large 
times, although the same nonlinearity can lead to large dynamic growth for early times. 



3 Properties of the Model Equation 

In this section, we will study some properties of the ID model equations. We will first 
consider the properties of some further simplified models obtained from these equations. 
Both numerical and analytical studies are presented for these simplified models. Based on 
the understanding of the simplified models, we prove the global existence of the inviscid 
Lagrangian model, which sheds useful light into our global existence analysis for the full ID 
model with or without viscosity. 



3.1 The ODE model 

To start with, we consider an ODE model by ignoring the convection and diffusion terms. 

(m)j = -2du (33) 
(v), = {uf - {v)\ (34) 

with initial condition 'S(O) = uq and 5(0) = vq. 

Clearly, if uq = 0, then u{t) = for all t > 0. In this case, the equation for v is decoupled 
from u completely, and will blow up in finite time if < 0. In fact, if vq < and uq is very 
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small, then the solution can experience very large growth dynamically. The growth can be 
made arbitrarily large if we choose uq to be arbitrarily small. However, the special nonlinear 
structure of the ODE system has an interesting cancellation property which has a stabilizing 
effect of the solution for large times. This is described by the following theorem. 

Theorem 2. Assume that Uq ^ 0. Then the solution {u{t),v{t)) of the ODE system 
(|33|) - (j34p exists for all times. Moreover, we have 

lim u{t) = 0, lim d{t) = 0. (35) 

>+oo t^+oo 



Proof. There are several ways to prove this theorem. The simplest way is to reformulate 
the problem in terms of complex variables^ . Let 

w = u + iv. 

Then the ODE system ()33|l - (j!^ is reduced to the following complex nonlinear ODE: 

dw 2 / \ 

— = iw , w{0) = Wo, (36) 



which can be solved analytically. The solution has the form 

1 — iwot 



In terms of the original variables, we have 



~U\ Uo{l + Vot) - UoVot 

«W = n , , .^2 ' (38) 



(1 + voty + {uoty 

vo{l + Vot) +uli 



v{t) = (39) 



It is clear from (jHHjl - ljH^ that the solution of the ODE system finn|) - (jMj) exists for all times 
and decays to zero as t — ^ +oo as long as Uq ^ 0. This completes the proof of Theorem 2. 



Remark 1. Note that the ODE model ()36|) has some similarity with the Constantin-Lax- 
Majda model which has the form ut = uH{u), where H is the Hilbert transform. By 
letting w = H{u) + iu and using the property of the Hilbert transform, Constantin-Lax- 
Majda show that their model can be written as the imaginary part of the complex ODE: 
Wt = \w'^- It is interesting to note that both models ignore the convection term and they 
have solutions that blow up at a finite time for initial condition satisfying ^(zo) = and 
H{u){zq) > for some Zq. However, as we will show later, the convection term plays an 
important role in stabilizing the ID model and should not be neglected in our study of the 

^We thank Prof. Tai-Ping Liu for suggesting the use of complex variables 



8 



Euler equations. By including the convection term in the ID model, we will show in section 
3.3 ans section 4 that no finite time blow-up can occur from smooth initial data. 

As we can see from fl38p - (jn^ . the solution can grow very fast in a very short time if Uq 
is small, but vo is large and negative. For example, if we let vq = — 1/e and = e for e > 
small, we obtain at t = e 

uie) = l/e^ v{e) = 1/e. 

We can see that within e time, u grows from its initial value of order e to 0(e^^), a factor of 
amplification. 

Remark 2. The key ingredient in obtaining the global existence in Theorem 2 is that the 
coefficient on the right hand side of is less than —1. For this ODE system, there are 
two distinguished phases. In the first phase, if v is negative and large in magnitude, but it is 
small, then v can experience tremendous dynamic growth, which is essentially governed by 

Vt = -v^. 

However, as v becomes very large and negative, it will induce a rapid growth in u. The 
nonlinear structure of the ODE system is such that u will eventually grow even faster than 
V and force (tt)^ — {vY < in the second phase. From this time on, v will increase in time 
and eventually become positive. Once v becomes positive, the nonlinear term, —v"^, becomes 
stabilizing for v. Similarly, the nonlinear term, —2uv, becomes stabilizing for u. This subtle 
dynamic stability property of the ODE system can be best illustrated by the phase diagram 
in Figure n 

In Appendix A, we prove the same result for a more general ODE system of the following 
form: 

(m)^ = -dvu (40) 
{v\ = {uf - {v)\ (41) 

for any constant d > 1. However, ii d < 1, it is possible to construct a family of solutions 
for the ODE systems (|i m) -PT |) which blow up in a finite time. 

3.2 The Reaction Diffusion Model 

In this subsection, we consider the reaction-diffusion system: 

(m)j = vu^^ - 2vu, (42) 
iv\ = vv,, + {uf - {v)\ (43) 

As we can see for the corresponding ODE system, the structure of the nonlinearity plays 
an essential role in obtaining global existence. Intuitively, one may think that the diffusion 
term would help to stabilize the dynamic growth induced by the nonlinear terms. However, 
because the nonlinear ODE system in the absence of viscosity is very unstable, the diffusion 
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Figure 1: The Phase Diagram for the ODE system. 



term can actually have a destabilizing effect. Below we will demonstrate this somewhat 
surprising fact through careful numerical experiments. 

In Figures ©-(jH), we plot a time sequence of solutions for the above reaction diffusion 
system with the following initial data 

Uo{z) = e(2 + sin(27r2;)), Vo{z) = sin(27r2;), 

e 

where e = 0.001. For this initial condition, the solution is periodic in z with period one. 
We use a pseudo-spectral method to discretize the coupled system ((121)- (03) in space and 
use the simple forward Euler discretization for the nonlinear terms and the backward Euler 
discretization for the diffusion term. In order to resolve the nearly singular solution structure, 
we use N = 32, 768 grid points with an adaptive time step satisfying 

Atn (I max{M"}| + I min{M"}| + | max{{5"}| + | min{5"}|) < 0.01, 

where m" and are the numerical solution at time t„ and t„ = tn~i + At„_i with the initial 
time stepsize Ato = 0.0 le. During the time iterations, the smallest time step is as small as 
0(10-1°). 

From Figure |21 we can see that the magnitude of the solution v increases rapidly by a 
factor of 150 within a very short time (t = 0.00099817). As the solution v becomes large and 
negative, the solution u increases much more rapidly than v. By time t = 0.0010042, u has 
increased to about 2.5 x 10® from its initial condition which is of magnitude 10~^. This is 
a factor of 2.5 x 10^^ increase. At this time, the minimum of w has reached —2 x 10®. Note 
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Figure 2: The solutions at t = 0.00099817, and t = 0.0010042, = 32768, z/ = 1. 



that since u has outgrown v in magnitude, the nonhnear term, u'^ — v"^, on the right hand 
side of the ^-equation has changed sign. This causes the solution v to split. By the time 
t = 0.001004314 (see Figure Ej), both u and v have split and settled down to two relatively 
stable traveling wave solutions. The wave on the left will travel to the left while the wave on 
the right will travel to the right. Due to the periodicity in z, the two traveling waves approach 
each other from the right side of the domain. The "collision" of these two traveling waves 
tends to annihilate each other. In particular, the negative part of v is effectively eliminated 
during this nonlinear interaction. By the time t = 0.00100603 (see Figure E)), the solution v 
becomes all positive. Once v becomes positive, the effect of nonlinearity becomes stabilizing 
for both u and v, as in the case of the ODE system. From then on, the solution decays 
rapidly. By t = 0.2007, the magnitude of u is as small as 5.2 x 10~^, and v becomes almost 
a constant function with value close to 5. From this time on, u is essentially decoupled from 
V and will decay like 0{l/t). 



3.3 The Lagrangian Convection Model 

Next, we consider the ID model equations in the absence of viscosity. The corresponding 
equations are given as follows: 

ut + '2'ipUz = -'2vu (44) 
Vt + 2^i)Vz = u^ -v'^ + c{t), (45) 

where v = and c(t) is defined in to ensure that vdz = 0. 
Introduce the Lagrangian flow map 

^^^ = 2^{z{a,t),t), (46) 
z{a,0) = a. (47) 
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Figure 3: The solutions at t = 0.001004314 and t = 0.001005862, N = 32768, z/ = 1. 




Figure 4: The solutions at t = 0.00100603 and t = 0.2007, N = 32768, u = 1. 
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Differentiating (^Ej) with respect to o;, we get 



dZfy^ dip 

— ^ = 2za — {z{a,t),t) = -2zav{z{a,t),t). 
at oz 

Denote f (a,t) = v{z{a^t)^t), u{a,t) = u{z{a,t),t), and J{a,t) = Za{(y,t). Then we can 
show that J, u, and v satisfy the following system of equations: 

^"^^"'^^ = -2J(«,t)t;(«,t), (48) 



dt 
du{a, t) 



-2u{a,t)v{a,t), (49) 



dt 

+ 3 I v^Jda- / u'Jda, (50) 
Jo 



dt 



with initial data J(a;, 0) = 1, u{a, 0) = Uo{a) and f (a, 0) = Voio:). Since f (2;, t)dz = 0, we 
have 

"1 

v{a,t)J{a,t)da = 0, (51) 





which implies that 

.1 .1 

J(Q;,t)(ia= / J{a,0)da = l. (52) 
io 

It is interesting to note that the ID model formulated in the Lagrangian coordinate retains 
some of the essential properties of the ODE system. In the following, we will explore the 
special nonlinear structure of the model equation to prove the global well-posedness of the 
ID model in the Lagrangian form. As we will see, the understanding of the ID model in the 
Lagrangian form gives critical insight in our understanding of the full ID model. 

Theorem 3. Assume that u{z,0) and v{z,0) are in C™[0, 1] with m > 1 and periodic with 
period 1. Then the solution {u, v) of the ID inviscid model will be in C™[0, 1] for all times. 

Proof. Differentiating the u and v equations with respect to a, we get 

dUn 



-2vUa — 2uVa, (53) 



dt 

= 2uua — 2vVa- (54) 
Multiplying by Ua and by Wq,, and adding the resulting equations, we have 

~{ul + vl) = -2v{ul + vl). (55) 

Therefore, we obtain 

^-^^\og{ul + vl)=-2v. (56) 
13 



Integrating from to t, we get 

= {VWa + ivo)l) e-'Io^^'^^s),s ^ (vT^^j^ + T«) ^(«,^), (57) 

where we have used 

which follows from PHj) and J{a,0) = 1. Using ()52p . we further obtain 

Vul + vl da < Wiuo)l + ivoya iU- / Jia,t)da=yiuo)l + {vo)l h^- (58) 



(59) 







In particular, we have 

/ \va\da< / ^/ul + vlda < \\^/{uo)l + {vo)l ||loo. 



^-1 ~ 

we get 



Since vJda = and J > 0, there exists ao(^) ^ [0, 1] such that v{aQ{t),t) = 0. Therefore, 



|{}(a,t)| = |/ iJ^da'\< \d{a',t)\da' < W + ivo)l iU- (60) 

This proves that 

Ml^ < Wi^o)l + {vo)l lUoo. (61) 

Using the equations for J and n, we also obtain 

e-2*^o < J(a,t) < e^*^", (62) 
< \\uq lU-e^*^", (63) 

where Cq = || VT^oH+TS^Ol 

The bound on J {a, t) in turn gives bound on Ua + v"^ through (jSTj). We can then bootstrap 
to obtain regularity of the solution in higher order norms. This completes the proof of 
Theorem 3. 

Next, we illustrate the behavior of the solution through numerical computations. We use a 
pseudo-spectral method to discretize in space and a second order Runge-Kutta discretization 
in time with an adaptive time-stepping. In Figures El and IHl we plot a sequence of snapshots of 
the solution for the inviscid model PH|) - ()50p in the Lagrangian coordinate using the following 
initial data 

u{a, 0) = 1, v{a, 0) = 1 - i exp-("-°•^)'/^ 



with e = 0.0001 and 6 = \fen . We can see that the solution experiences a similar splitting 
process as in the reaction diffusion model. In Figure [71 we perform a similar computation in 
the Eulerian coordinate with e = 0.00001. We can see that as the solution v grows large and 
negative, the initial sharp profile of v becomes wider and smoother. This is a consequence 
of the incompressibility of the fluid flow. If we change the sign of the convection velocity 
from 2^ to —2^, the profile of v becomes focused dynamically and develops an unphysical 
"shock-like" solution, which seems to evolve into a finite time blowup, see Figure |H1 
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Initial condition for v, N=4096 



500 1000 1500 2000 2500 3000 3500 4000 4500 



-150 - 

-200 - 
-250 - 



Solution V at t=0.01 88, N=4096, viscosity=0 



500 1000 1500 2000 2500 3000 3500 4000 4500 



Figure 5: The Lagrangian solution at t = and t = 0.0188, = 4096. 



Solution V at t=0.023, N=4096, viscosity = 



Solution V at t=0.0337, N=4096, viscosity = 




1600 1700 1800 



2000 21 00 2200 2300 2400 2500 2600 



1600 1700 1800 



2000 21 00 2200 2300 2400 2500 2600 



Figure 6: The Lagrangian solution at t = 0.023 and t = 0.0337, N = 4096. 
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Figure 7: The Eulerian solution, N = 4096. 



Time sequence of v from t=0 to 0.0033, N=8192, wrong convection sign, vis=0 




-450 - 

-500' ' ' ' ' ' 

3950 4000 4050 4100 4150 4200 4250 



Figure 8: The Lagrangian solutions with the wrong sign, = 4096. 
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4 Global Well-Posedness of the ID Viscous Model 



Based on the understanding we have gained from the previous sections, we are ready to 
present a complete proof of the global well-posedness of the full ID model. It is not easy to 
obtain global regularity of the ID model by using an energy type of estimates. If we multiply 
the M-equation by u, and the -y-equation by v, and integrate over z, we would arrive at 

-— / u^dz = — 3 / {uYvdz — y I u^dz, (64) 

1 d 
2dt 



JO JO 

1 pI pi pi 

v'^dz = I u^vdz — 3 / {v^dz — v^dz. (65) 
Jo Jo Jo 



Even for this ID model, the energy estimate shares the some essential difficulty as the 
3D Navier-Stokes equations. It is not clear how to control the nonlinear vortex stretching 
like terms by the diffusion terms. On the other hand, if we assume that 

\\v\\Loc,dt < OO, 





similar to the Beale-Kato-Majda non-blowup condition for vorticity then one can easily 
show that there is no blow-up up to t = T. 

In order to obtain the global regularity of the ID model, we need to use a local estimate. 
We will prove that if the initial conditions for u and v are in with m > 1, then the 
solution will remain in for all times. 

Theorem 4. Assume that u{z,0) and v{z,0) are in C"^[0, 1] with m > 1 and periodic with 
period 1. Then the solution {u, v) of the ID model will be in C""[0, 1] for all times. 

Proof. Motivated by our analysis for the inviscid Lagrangian model, we will try to obtain 
a priori estimate for the nonlinear term + v"^. Differentiating the w-equation and the 
S-equation with respect to z, we get 

{u^)t + 2i){u^)^ - 2vu^ = -2u^v - 2uv^ + //(m^)^^, (66) 
{v^)t + 2%i){vz)z - 2vvz = 2uuz - 2vv^ + i'{vz)zz- (67) 

Note that one of the nonlinear terms resulting from differentiating the convection term 
cancels one of the nonlinear terms on the right hand side. After canceling the same nonlinear 
term from both sides, we obtain 

{u^)t + 2i}{u^)^ = -2uvz + i'{u^)zz, (68) 

{V:,)t + 2i){V:,)^ = 2UU^ + v{V:,)^^. (69) 

Multiplying by Uz and (jU^ by -u^, we have 

]^{ul)t + ij{ul)z = -2uUzVz + uUz{iMz)zz, (70) 
^{vl)t + i'{vl)z = 2uuzVz + yvz{vz)zz- (71) 
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Now, we add (f7n|) to ()7H). Surprisingly, the nonlinear vortex stretching-like terms cancel 
each other. We get 

{^l + vl)^ + 2i> {ul + v])^ = 2u [u,[u,),, + v,(v,),,) . (72) 
Further, we note that 

(m^ + f;^)^^ = {2U^U^^ + 2v^Vzz)^ = 2 {Uz{Uz)zz + Vz{Vz)zz) + 2 [{Uzzf + [Vzzf] ■ 

Therefore, equation ()72|1 can be rewritten as 

{ul + v^;)^ + 2V^ {ul + v^^)^ = y{ul+ il)^^ - 2v [(w,,)^ + {vz-:f\ . (73) 

Thus, the nonlinear quantity, (u^ + ^z)? satisfies a maximum principle which holds for 
both z/ = and z/ > 0: 

WuI^vIWl^ < ||(mo)' + (^^o)'||l-. (T4) 
Since v has zero mean, the Poincare inequality implies that < Co, with Co defined by 

Co =11 {{uo)l + iVo)iy'\\L^. 

The boundedness of u follows from the bound on v: \\u(t)\\L°° < HmoHl"" exp(2Coi). The 
higher order regularity follows from the standard estimates. This proves Theorem 4. 



5 Construction of a family of globally smooth solutions 

In this section, we will use the solution from the ID model to construct a family of globally 
smooth solutions for the 3D axisymmetric Navier-Stokes equations with smooth initial data 
of finite energy. We remark that a special feature of this family of globally smooth solutions 
is that the solution can potentially develop very large dynamic growth and it violates the 
so-called smallness condition required by classical global existence results jHllZIl- 

Let Uilzjt), u)i{z,t), and ip be the solution of the ID model problem. We will construct 
a family of globally smooth solutions of the 3D Navier-Stokes equations from the solution 
of the ID model problem. Denote by u{r,z,t), Co{r,z,t) and ip{r,z,t) the solution of the 
corresponding 3D Navier-Stokes equations. Further, we define 

Ul = u/r, Ul = ijj/r, ipi = ip/r. (75) 

Let 0(r) = (pQ^r/Ro) be a smooth cut-off function, where 0o(r) satisfies 0o(r) = 1 if < r < 
1/2 and and 0o(^) = if r > 1. Our strategy is to construct a family of globally smooth 
functions Ui, ui and ipi, which are periodic in z, such that 

u = r {ui{z,t)(f){r) + ui{r, z,t)) = u + u, (76) 
u = r{uji{z,t)(l){r)+uji{r,z,t)) =1}) + ^), (77) 
^ = r{ilJi{z,t)(j){r) + xlji{r,z,t)) =00 + 00, (78) 
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is a solution of the 3D Navier-Stokes equations. 

With the above definition, we can deduce the other two velocity components v'^ and 
as follows: 

= —ipz = —f^(p{''^)4'iz + v'^{r, z, t) = v"^ + , (79) 

v' = = 0(2^/'i) + r(j)ripi + v'{r, z, t) = v' + v'\ (80) 

r 

With the above notations, we can write the velocity vector into two pars as u = u + u. 
We will choose the initial data for the ID model of the following form: 

^i(^'0) = «i(^,0) = ^f/i(zM), u:,{z,Q) = AW,{zM), (81) 

where A and M are some positive constants, '^i{y), and Ui{y) are smooth periodic functions 
in y with period 1. Moreover, we assume that and Ui are odd functions in y. Clearly we 
have Wi = —{'^i)yy, which is also a smooth, periodic, and odd function in y. It is easy to see 
that this feature of the initial data is preserved by the solution dynamically. In particular, 
ipi{z,t), Ui{z,t), ui{z,t) are periodic functions in z with period 1/M and odd in z within 
each period. Using this property and the a priori estimate ()74|1 . we obtain the following 
estimate for the solution of the ID model: 

\\Mt)\\L^<Co^, (82) 

\\u,mL- < Co^, UUt)\\L- < Co^, (83) 

||cui(t)|U=o < CoA, \\u,M\l^ < CoA, (84) 

where 

Co= II {Ul + W^y^U^. (85) 

Remark 3. As we know from the discussions in the previous sections and as indicated 
by (jH^ - IH3j) . if the regularity of the periodic profiles in the initial condition, i.e., Ui and 
\I'i, is very poor, the solution Ui{z,t) and il)z{z,t) will grow very fast dynamically. The 
amplification factor is determined by Co defined in (jH3j) . 

Let Ro = M3. From (jH^ - (jS3|) and the definition of u, we have 

\\u\\l2^ ARl/M = A/^, \\V\y\\L2^ ARl = A^. (86) 

We would like to emphasize that the corresponding 3D solution defined by (ffH|) - (ffH|) 
in general does not preserve the same special structure in the z direction of the ID model 
problem since the correction terms, Ui, uji and ipi, are periodic in z with period 1 instead of 
period 1/M. 
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We assume that the initial conditions for mi, tui, and ipi are chosen in such a way that the 
principal contributions to the energy and the enstrophy come from u, the mollified solution 
of the ID model. Specifically, we assume that the initial condition for u satisfies: 



\\uo\\ ^ AR'q/M = A/ VM, ||Vuo||l2 ~ = ^vM. (87) 

Thus, we have 

I|uo||l2||VCio||l2 ~ A'. (88) 

By choosing A large enough, the above product can be made arbitrarily large. Thus it 
violates the classical "smallness" condition that guarantees the global existence of the 3D 
Navier-Stokes equations PT] . 

Furthermore, we have from the energy inequality that 

Mt)\\L2 <\\uo\\l^ <A/Vm. (89) 

Using the above bound and ()86|). we obtain a priori bound for the perturbed velocity field, 
u in norm: 

||u(t)|U. <^. (90) 

Let f = u1, and define 

H^{t) = j {f + ujl)rdrdz = j {u\ + ujl)rdrdz, (91) 
E^{t) = j {\\/f\^ +\VuJi\^)rdrdz, (92) 

where the integration is over [0, 1] x [0, oo). 

If we further assume that the initial conditions for ui, ui, and ipi are odd functions of z, 
then it is easy to verify that ui, ui and ipi are odd functions of z for all times. Since ui, uji, 
and ipi are also odd functions of z, we conclude that Ui, ui, and ipi are odd functions of z 
for all times. It follows by the Poincare inequality that we have 



frdrdz < I flrdrdz < I \Vf\^rdrdz, (93) 



This implies that 



H <E. (95) 
Now we can state the main theorem of this section. 

Theorem 5. Assume that the initial conditions for mi, uji and ipi are smooth functions 
of compact support and odd in z. For any given A > 1, Co > 1 and z/ > 0, there exists 
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C{A,Co,u) > such that if M > C(A,Co,z/) and H{0) < 1, then the solution of the 3D 
Navier-Stokes equations given by (f7Bj) -()78 p remains smooth for all times. 

Proof. First of all, we can use p3|) - (fT3j) to derive the corresponding evolution equations for 
Ml, uJi and ipi as follows: 



J. 

(a)l)i + v''{uJi)r + V^'i^UJi)-, = (ul)^ + V \ + Lb^rr + , (97) 

Substituting (fTUj) into and using (fTSj) . we obtain an evolution equation for Ui. 

^ h 'u'''Uir + "C^Mlz = Z/Ami + 2'?/'i2'Ui — U\i^ — v'^UiCpr (99) 

- 0-u^-ui2 + Z^A(-Ui0), 
where we have used A to denote the modified Laplacian operator defined by 

Aw = Wzz + Wrr H = W^z + ArW. 

r 

On the other hand, we know that ui satisfies the ID model equation: 

Ult + 2V^iMi^ = z/Mi^^ + 24}uUi. (100) 
Multiplying ()100|1 by and subtracting the resulting equation from (j99|) . we have 

Uu + v'^Uir + v'^Uiz = Z/Aui + 2 (^V^lz'Si - (p^jjlzUi^ (101) 

— V^'Ui(j)r — (['"^r + 2(0 — 1)]'01 + f ^) + Z/'UiA^0. 

Similarly, we obtain 

cjif + ■u'^cuir + -y^cui^ = uAuJi + [{ui + u(t))l ~ ul^(f)) (102) 

— ?)''ti}i0r. — ([r0r + 2(0 — 1)]'01 + f^) tUlz + ULJiArCf). 

We divide the analysis into two parts. The first part is devoted to estimates of the 
velocity equation. The second part is devoted to estimates of the vorticity equation. 

Part I. Estimates for the velocity equation. 

First we will present our analysis for the velocity equation. 
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Multiply pOlj) by uf and integrate over [0, 1] x [0, oo). Using the incompressibility con- 
dition 



we get 

Id f A J J 3z/ f 1 , 2\|2 

/ u^rdraz < / V Mi raraz 

Adtj ^ - 4 y I ^ 

+ 2U^..^.-H..u;)rdrd.-Uu.^.U^drd. 

hasthesameorder = — — / |V('Ui)| rdrdz + I + II + III + IV, (103) 



where we have used the fact that 



ulAuirdrdz = I t^i { "^izz + - — — ^ H — ) rdrdz 



((^i)^)^ + ((«i)r)^ rdrdz / u\{S\,z,t)dz 



3 
~2 
_3 
~4 

- "i / [(^^')^)' + ^^^^^ • (1°^) 

In the following, we will estimate the right hand side of ()1()3|1 term by term. 

Estimate for the I-term. 
Using (UnD-CZHD, we have 

/ = 2 y Ml {ipizUi + + (t)Uii)u + (0^ - rdrdz 

= la + h + Ic + Id. (105) 

Using the Holder inequality, we have 

Ia<2\MLA\f\\l^- 

Note that 

llV^i^ll^a = J i)l^rdrdz = J i)l^d{r'^ /2)dz = ~^ j ripuipurrdrdz 
1 A 

< -Wf^^lzhAl'^lzrh^ < rll^lzrlUa, (106) 

Z Ma 
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where we have used ripiz = V'z cind (plj) to obtain 

A 



Ms' 



On the other hand, using the Sobolev interpolation inequahty, we have 

ii/iiL.<ii/iiiiiv/iii. 

This imphes that 

Mi M- 

where C2 = ||Ar0o||L°°5 and we have used 

II , II ^ C2C0A 

Mi 

which we prove in Appendix B. 
The estimate for follows from 

h < 



2Co— [ u^rdrdz < 2Co—H\ 



As for Jc, we use (jHH|) . (jl(J6|) . and the Holder inequality to obtain 

ic < 2Co^\ML4f\\l 

< 2Co^^||^i..||i||/||i||V/||i 

3 ^ 



'2^2 

Mi Mi 



2Co(A)2 5 3 2V^CoM^,3_3 



where we have used ()1()7|1 and the Sobolev interpolation inequality 

ii/iiL3<coii/iiiiiv/iii. 

Finally, we use (jH^ and the Holder inequality that 



Therefore, we obtain 



/ < 2^HEi + 2^^^^HlEi + 2cAH^ 
Mi Mir M 



3 3 



Mi Mf Af2-i/8 
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Estimate for the Il-term. Using (jH^ . and tlie Holder inequality, we have 



CiCqA 3i 7 7 ^ CiCqA 



3 



< ''^^yWhW'^fWhi'^^^HiEi, (112) 



^^;:.';iliiv/iiir-2:^^ 



where ci = || (0o)r|U°°, and we have used the Sobolev interpolation inequality ()109|) . 



Estimate for the Ill-term. Using ()82|) . ()84|) . and following the same steps as in our 
estimate for the J^-term and the //-term, we get 

HI < + / \u\\rdrdz + CoA / \v'\\u\\rdrdz 

Jr<Ro J 



^ i^ + ^i)T;ihmH^ + CoA\\v^h4f\\h 



'jVf 2-1/8 

C^A^ ,,3 , CoA'' 



n3 „3 



£ P + '^OjprTTi^f + (113) 



Estimate for the IV-term. 

Using (jH!^ and the Holder inequality, we have 



VC2CqA f I i3 J J ^ VC2CqA ^ VC2C0A ^^3 

W 



IV < I \uA^rdrdz < —^-^- fft ps < ni4) 



Part II. Estimates for the vorticity equation. 

Next, we will present our analysis for the vorticity equation. Multiplying (jlU2|) by uji 
and integrating over [0, 1] x [0, 00), we get 

(p [[rcpr + 2(0 — + f ^) uJizUJirdrdz + cuiArCpuJirdrdz 



V \ \VuJiYrdrdz + I + 11 + III + IV, (115) 



where A^^ = 0T-r + —■ 



We will estimate the terms I to JV^ one by one. 
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Estimate for the /-term. 

Using (fTUj) and (jHHj) and integration by parts, we have 

/ = — J {ul — u\(j))uoizrdrdz = — j {u\ + 2ui(j)Ui + (0^ — (f))u\)uoizrdrdz 

< / u^\uji^\rdrdz ^ — — / \ui\\uJiz\rdrdz ^ — — ||u;i^||l2 

< ((/ uir,r,.) ' + 2^ ' + ] M.. (116) 

Let r = ru. It is easy to show that F satisfies the following evolution equation (see also 

EDI) 

r 

Moreover, for u smooth, we have T\r=o = 0. Thus, F has a maximum principle, i.e. 

l|r||L°° ^ l|rol|L°° ^ cq. 

This implies that 

|2 |^|~|i 2|j-|^ id2^0^ . Co^~ 

IVI ]\/[ 2 

Therefore, we obtain 

u\rdrdz = I {ul)'^d{r'^ /2)dz = — rul{ul)rrdrdz 



< r u^rdrdz j ||V/||l2 

< \\rui\\l2 (^j r^lrdrdz^ ' ||V/||l2 

< So^(/ ulrdrdzyWVfh^. 



On the other hand, we have 



J ulrdrdzy < y utrdrdzy WfWl < (^j u\rdrd^' \\f\\^,\\W f\\f,. 



Combining the above estimates, we obtain 



1 ~- 1 

2 Cn A 7 3 _ 19 



u^.rdrdz] <^^H^E^. (117) 
'Mr 
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Thus, we obtain 



1 



I < ^^H^E^ + / u^rdrdz + 

' Mr M ^ \J ^ J M- 



4 



Mr Ms Mi 



Estimate for the JJ-term. 

Using (jH^ . (jnOJ, and the Holder inequahty, we get 

— ciCoA 

-fin 



^0 M2 M4 



Estimate for the ///-term. 

Integration by parts gives 

777 = j (j){[r(f)r + 2{4) - + v^^)iDiUJirdrdz 

+ j (t){[r(l)r + 2{(t)-l)]i)i+v'')uJiUJi^rdrdz 
We first study the term J (pvluiUirdrdz. Note that using (jH^ . we have 

On the other hand, we have by the Sobolev interpolation inequality that 



h\\V(<t>v')A\h 
< ^l|V(r«),||i 

Mi 

which we prove in Appendix B. 



where we have used 
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Thus we use (jHSl) that 



III < {2 + ci)^I^\\ui\\l2+ [ ,\,^^ HE-^+ ^ I H 



+ (2 + ci)^^i?o||o;i.|U2+CoA|| 



/1 2 or'n A i 1 

< (2 + ci + 2^2)^H + ,° ^ + (3 + c] 



M3 



(^2/12 
M2 



(121) 



Estimate for the IV-term. 

Let g{z,t) = LUi{r],t)dr]. Then we have = uJi and < Thus we have 

= u J Cui{Ar(f))uJirdrdz ~ ^ J gzi^r<P)^irdrdz 

= -u [ g{Ar(l))uJizrdrdz < [ \uJi^\rdrdz 

J MKq Jr<R,, 



UC2C0A „ ,, ,, ^ UC2C0A 



By adding the estimates for J u\rdrdz to those for J ufrdrdz, we obtain an estimate for 
^H^- Note that except for the diffusion terms, each term in our estimates from I to IV can 
be bounded by 

where g{II) is a polynomial of if with positive rational exponents and positive coefficients 
that depend on Co, A, and z/, and e = for some 7 > 0. Putting all the estimates together, 
we get 

jH' < -'^E^ + eg{H) < -'^H^ + egiH), (123) 

since H < E. 

For given A > 1, Cq > 1, and z/ > 0, we can choose M large enough so that 

-^ + e^7(l)<0. 

Thus, if the initial condition for ui, uoi and ipi are chosen such that -ff(O) < 1, then we must 
have 

H{t) < 1, for all t > 0. 

Using this apriori estimate on II{t), we can easily follow the standard argument to prove 
the global regularity of ipi, Ui and uji in higher order norms. This completes the proof of 
Theorem 5. 
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Appendix A. 



In this appendix, we prove the following result for the generalized ODE system. 



Theorem A. Assume that Mq 7^ and d > 1. Then the solution {u{t),v{t)) of the ODE 
system f)40|) - (PT|) exists for all times. Moreover, we have 

lim u{t) = 0, lim v{t) = 0. (124) 



Proof. We first make a change of variables into the polar coordinate ^ 

V = r cos 6, u = r sin 6. (125) 

Substituting the above change of variables into the ODE system, we obtain 

r' cos e - r(sin 6)6' = sin^ 6 - cos^ 6, (126) 
r' sin 6 + r(cos 6)6' = —dr^ cos 6 sin 6. (127) 

From the above equations, we can easily derive 

r' = -r^ cos 6 (cos^ 6 + {d-l) sin^ 6) , (128) 
6' = -rsm6 {{d-l)cos^ 6 + 8111^6) . (129) 

Note that if Uq > 0, then u{t) > as long as | jQv{s)ds\ < 00. Similarly, if uq < 0, then 
u{t) < 0. Thus, if the solution starts from the upper (or lower) half plane, it will stay in 
the upper (or lower) half plane. Without loss of generality, we may consider the solution 
starting from the upper half plane. It follows from p29|l that 6' < since d > 1, and r > 0. 
Therefore, 6{t) is monotonically decreasing. On the other hand, 6{t) is bounded from below 
by zero. As a result, the limit of 6{t) as t —>■ 00 exists. Let us denote the limiting value as 
6. Clearly, we must have 

lim 6' = 0, lim 6 = 6. (130) 

First, we consider the case that the solution starts from the second quarter (y axis 
included). We claim that this solution must cross the y-axis into the first quarter. If the 
solution stay in the second quarter forever, then 6 must be no less than 7r/2. From p29p 
and ()13U|) . we know that 

lim r = 0. (131) 

However, from ()128|) . we have r' > 0, which contradicts with p3H) . The contradiction 
implies that the solution must cross the y-axis at a later time. 

Now we only need to consider the case when the solution starts from the first quarter 
since the system is autonomous. Since 6{t) decreases monotonically, we obtain 

r' < -cos%6o)r^, (132) 
^This proof was inspired by a discussion with Mr. Mulin Cheng. 
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where we have used the fact that d> 1 and cos^ 9 + (d — l) sin^ 9 > cos^ 9. Solving the above 
ODE inequahty gives 

^it) < 1 T^^TT^- (133) 

^ ^ - l + ro(cos3^o)i 

Thus, we conclude that 

lim r{t) = 0. (134) 

To determine that limiting angle, 9, we use the fact that 

'u'\ dtan{9) 



= it) 



tan^{9) - 1 

Since d > 1, we conclude that ^ = 0, which implies 



lim 9{t) = 0. (135) 
This completes the proof of Theorem A. 

Appendix B. 

In this appendix, we prove the following two estimates which relate the norm of the 
derivatives of ip to that of uJi. 

Wi^lzzh'^ + llV^lrzlUa + ll^irrlUa + || — ||l2 < ||wi||l2 + (136) 

r Mi 

and 

WvLWl^ = l|V(^ + ^i..)|U2 < \\Vw,h2 + ^^, (137) 
r Mi 

where C2 = || Ar</)o||L°=. 

Proof. From the definition, we have 

—Atpi = wi. 

Using the definition of ipi and wi, we can rewrite the above equation as 

-wi = A^i + (A,0)V5i. (138) 
Multiplying (jl38|) by tpi^z and integrating over [0, 1] x [0, oo), we obtain: 

IkilUHlV^iz^lU^ > J (AipiiJizz - {Ar(l))i}iiJi;,;,)rdrdz (139) 

> J ii'izz + ^irz^drdz - 2 J ipirzi'izdrdz - J \^izz\rdrdz (140) 



> I irizz + rirz)rdrdz+ I ^|Jl(0,z,t)dz-J^J\^P^,,U2, (141) 
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where we have used (|S^. This imphes that 



C2C A 

11^12^11^2 + ||^lr-2||L2 < ||Wi||l2 + ^ ° ■ (142) 

Mi 

Next, we muhiply p38|l by /S.ripi and integrate over [0, 1] x [0, 00). We obtain by using 
a similar argument that 

||wi||i2||A^^i||i2 > j (A^iA^V^i - {/^rffyi'i^Mrdrdz (143) 



> / [(A,.V^i)' + V^L]rrfrrfz- ^||A.V^i|U2. (144) 



On the other hand, we note that 



{Arijjiyrdrdz = I (ipl^.^ + 9'^^)rdrdz + Q j ipir^irrdrdz 



J {ipfrr + '^^^ydrdz ^ 3 y '?Air(0, z, t)dz 



where we have used the fact that ipiriO, z,t) = since ipir is odd in r. Thus we obtain 

||^lrr||L2 + || — ||l2 < ||A^V'i|| < \\Wi\\l2 + (145) 

r Mi 

Combining estimate p44|) with p45|) gives the desired estimate ()136p . Similarly, we can 
prove p37|) . 
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